Nanotubes in polar media: polarons and excitons on a cylinder 
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Electrons and holes on semiconducting nanotubes immersed in sluggish polar media can undergo 
self-localization into polaronic states. We evaluate the binding energy ol of adiabatic Frohlich- 
Pekar polarons confined to a cylindrical surface and compare it to the corresponding exciton binding 
energy _E£ XC . The ratio E^° l /El xc is found to be a non-monotonic function of the cylinder radius R 
that can reach values of about 0.35, substantially larger than values of about 0.2 for 2d or 3d systems. 
We argue that these findings represent a more general crossover effect that could manifest itself in 
other semiconductor nanostructures in 3d polar environments. As a result of the strong polaronic 
effect, the activation energy of exciton dissociation into polaron pairs is significantly reduced which 
may lead to enhanced charge separation. 
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PACS numbers: 73.22.-f, 73.21.-b, 71.35.-y, 71.38.-k 



I. INTRODUCTION 



Low-dimensional semiconductor structures such as 
quantum wells, quantum wires, nanotubes and conju- 
gated polymers are important for practical applications 
and interesting scientifically. It is known that the con- 
finement of the motion of charge carriers in some direc- 
tions leads to increased effects of the Coulomb interaction 
on system excitations. In this paper we are concerned 
with two types of such effects: excitonic and polaronic. 
The excitonic effect refers to the formation of Coulombi- 
cally bound electron-hole pairs, Wannier-Mott excitons, 
which progressively affect optical properties of semicon- 
ductors: the exciton binding energy E^ xc increases from 
its 3d value to 2d and, further on, to quasi-ld magnitudes 
(E£ xc diverges in pure ld)i* The polaronic effect occurs 
in polar media, where the Coulomb field of an individual 
charge carrier causes the polarization (deformation) of 
its surroundings resulting in the carrier self-localization 
into polaronic states. Such polarons have been especially 
extensively studied in the context of 3d ionic crystals and 
polar semiconductorsi*ii** The polaronic effect also grows 
with the confinement: the polaron binding energy E^ Ql 
in 2d is larger than in 3d and would diverge in pure Id** 
Various aspects of the excitonic and polaronic effects have 
been explored for many specific low-dimensional systems. 



Of particular interest to us is a relationship between 
and E%° 1 , each of the binding energies understood 
as measured from the band edges in the absence of the 
polaronic effect. The ratio E^° l /E£ xc of the binding en- 
ergies has a clear significance for the relative energetics 
of closely-bound and well-separated electron-hole pairs 
that is expected to affect practically important processes 
of charge separation and recombination. The bare value 
j^exc signifies the ionization energy, whether thermal or 
photo, required for "unbinding" of the exciton into well- 
separated band-edge electron and hole. In the presence 
of the polaronic effect, however, the thermal ionization 



(dissociation) would occur into distant electron-polaron 
and hole-polaron so that the exciton thermal ionization 
energy is reduced from E£ xc to E^ xc —2E^° l . We are inter- 
ested in how much of this relative reduction might be pos- 
sible to achieve due to the formation of strong- coupling 
(adiabatic) polarons. Our discussion here is restricted 
to systems with equal electron and hole effective masses 
so that the electron-polaron and hole-polaron have the 
same binding energies while the neutral exciton does not 
experience the adiabatic polaronic energy correction^* 

It is instructive to look at the ratio E^ ol /E^ xc based on 
the results known for isotropic systems of "well-defined" 
dimensionality, that is, for purely 3d, 2d and Id systems. 
One would then find that, while each of the binding en- 
ergies increases with more confinement, their growth oc- 
curs nearly "in proportion" so that the ratio does not 
change significantly. Indeed, classic Pekar's result* for 
3d adiabatic polarons would translate into the maximum 
ratio of about 0.22. The exciton binding in 2d increases 
by a factor of 4 from its 3d value* but almost so does 
the polaron binding energy in 2d* resulting in the ratio 
~ 0.20. Moreover, if the divergent purely Id binding en- 
ergies are taken (parametrically) for their ratio, then the 
result of Ref. 6 would translate into maximum E^ ol /E^ xc 
of approximately 0.17. 

From the standpoint of this data, the results of our 
recent quasi-ld calculations** for polarons and excitons 
on nanotubes immersed in a 3d polar medium yielding 
Efr 0l /E% xc in excess of 0.3 appear quite surprising. In 
the particular case of the tubular geometry, charge car- 
riers are confined to the motion on a cylindrical surface. 
Our study— was restricted to relatively small cylinder 
radii R. In this paper we will use a direct variational 
approach to calculate adiabatic polarons for arbitrary 
R thereby enabling an assessment of the evolution of 
the ratio El° l /El xc between the purely 2d (R -» oo) 
and quasi-ld regimes. The corresponding calculations 
for excitons on a cylindrical surface have been performed 
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recently^ii^ which are very much in line with our exci- 
ton data. We will explicitly demonstrate (see Fig.^b)) 
that the ratio of the polaron and exciton binding energies 
exhibits a non-monotonic behavior as a function of the 
cylinder radius and can achieve as large values as about 
0.35 at intermediate R where the cylinder circumference 
is roughly comparable to an appropriate Bohr radius. We 
believe that our demonstration of E?° /E% xc ratios above 
the values in purely 3d and 2c? systems can be rationalized 
by invoking simple physical arguments. These arguments 
also suggest that the found relative increase of the pola- 
ronic effect may be reflective of a general behavior that 
might be characterized as a crossover effect. 

Consider a gradual increase of the confinement, e. g., 
by starting to decrease the radius of a very large cylin- 
der in going from the purely 2d system towards quasi-ld 
or by starting to decrease the thickness of a very thick 
quantum well in going from the purely 3d system to- 
wards 2d. Simplistically stated, the spatial size (extent) 
of a yet unconfined polaron wave function is larger than 
that of an unconfined exciton. As the confinement in- 
creases, therefore, the polaron can start experiencing a 
substantial growth of its binding energy due to the con- 
finement before a "proportionally" substantial growth of 
the exciton binding energy sets in. With still further in- 
crease of the confinement, both polaron and exciton bind- 
ing energies will be reflecting fuller confinement effects 
resulting in the corresponding trend of the decreasing 
ratio E^ ol /E^ xc . Of course, quantitative aspects of the 
evolution can vary for different systems and need to be 
evaluated accordingly. We also note that the above con- 
sideration tacitly assumed the existence of a uniform 3d 
polarizable medium with the confinement affecting only 
the motion of charge carriers. Strong violations of this 
assumption can significantly affect the outcome for the 
ratio of binding energies. 

While serving as a suggestive illustration of a possibly 
general behavior, it is polarons and the relationship of po- 
larons and excitons on a cylindrical surface that are the 
subjects of our direct interest in this paper. Organic and 
inorganic tubular (nano)structures attract a great deal of 
attention and are considered candidate systems for im- 
portant applications like (photo)electrochemical energy 
conversion and (photo)catalytic production and storage 
of hydrogen as well as in optoelectronics. On one hand, 
their extended size along the tube axis can facilitate 
very good electron transport in that direction. On the 
other hand, tubes can expose large areas of both exte- 
rior and interior surfaces to facilitate surface-dependent 
reactions. Many of these applications involve contact 
with polar liquid environments such as common aqueous 
and non-aqueous solvents and electrolytic solutions which 
can provide conditions appropriate for the strong pola- 
ronic effect^ - thereby changing the nature of charge car- 
riers. A widely known example of the tubular geometry 
is single- walled carbon nanotubes (SWNTs) and, in fact, 
"chiral selective reactivity and redox chemistry of carbon 
nanotubes are two emerging fields of nanoscience" i& We 



note that the importance of the excitonic effects in the 
optics of semiconducting SWNTs is well established now 
with the radius independent binding energies E^ xc exper- 
imentally measured in some SWNTs to be in the range 
of 0.4 — 0.6 eVi^i^i^ There is also a growing evidence 
of the environmental effects on the electronic properties 
of SWNTs JLiUiL The assumption of equal electron and 
hole effective masses is a good approximation for SWNTs. 

It can be anticipated that the polaronic effect would 
have an influence both on charge-transfer reactions and 
charge carrier dynamics on the tubes. In the context 
of our discussion of the relative energetics, a substantial 
reduction of the activation energy due to the polaronic 
effect should be expected for electric-field-assisted exci- 
ton dissociation and charge separation on semiconducting 
nanotubes. 



II. EXCITON AND POLARON ENERGY 
FUNCTION ALS 

Following Refs. OandlH the basic model we consider 
in this paper assumes that electron and hole are particles 
of the same effective mass m e = rrih = m whose motion is 
restricted to the surface of the cylinder of radius R. The 
cylindrical surface itself is immersed in the uniform 3d di- 
electric continuum characterized, as is common in studies 
of the polaron2i&4 and solvation 20 effects, by two mag- 
nitudes of the dielectric permittivity: the high-frequency 
(optical) value of eoo and the low-frequency (static) value 
of e s . In the case of liquid polar media, the slow compo- 
nent of the polarization is ordinarily associated with the 
orientational polarization of the solvent dipoles and it is 
typical 10 20 that e s 3> £oo- The fast component of the 
polarization follows charge carriers instantaneously. The 
slow component of the polarization, on the other hand, 
is considered static in the adiabatic picture to determine 
the electronic states; the slow component then responds 
to the averaged electronic charge distribution. 

With m e = m/j, there is no net charge density asso- 
ciated with the ground state of the neutral exciton that 
would cause a static polarization of the slow component 
("non-polarizing exciton"). The Coulomb interaction be- 
tween the electron and the hole in the exciton is screened 
only by the high-frequency dielectric response eooi&Si 
The total energy of the exciton as a functional of the 
normalized wave function ip{v) of the electron-hole rela- 
tive motion (reduced mass to/2) is then 

£ c * c {iP} = x cxc -J7 0XC 

whose global minimum determines the exciton ground 
state wave function and its binding energy. For our ge- 
ometry, position vector r = (x, y) is confined to the cylin- 
drical surface, where we choose x to be along the cylin- 
der axis and — ttR < y < nR along the circumferential 
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direction; \Vip(r)\ 2 = (dtp/dx) 2 + (dtp/dy) 2 . Coulomb 
interaction is determined by the physical distance in the 
3d space; in the flat geometry it would be D(r) = |r|, for 
the points on the cylindrical surface 



D(r) = 



1/2 



(2) 



Optimization of the functional for the ground state of 
the exciton on a cylinder has been performed in Refs. ^3 
and ^2 with the results in a very good agreement with 
our data to be used in the comparison with the polaron. 

Polarons we discuss here are of the large-radius 
Frohlich-Pekar type where the Coulomb field of an in- 
dividual charge carrier supports a self-consistent dielec- 
tric polarization pattern surrounding the carrier. The 
formation of large-radius polarons (self-localization, self- 
trapping) occurs due to the interaction of charge carriers 
with the "slow" component of polarization. As the fast 
component of polarization does not contribute to the po- 
laronic effect, it is the effective dielectric constant e*: 



1/e* = 1/eoo - l/e s , 



(3) 



that affects the coupling strengthi 2 *^ In the adiabatic 
approximation, the normalized wave functions ip( r ) °f a 
self-localized charge carrier correspond to the minima of 
the following polaron energy functional: 

E pol {ip} = K pol -U po1 



|- /dr |VV(r)| 2 
2m 



e 

2e* 



dri dr 2 



\^{vi)\ 2 |^(r 2 )| 2 
I>(ri-r 2 ) 



(4) 



whose global minimum we will be seeking for the ground 
state of the polaron. The functional I0J is a result of the 
optimization of the total adiabatic energy functional with 
respect to the polarization of the medium thereby ex- 
hibiting an effective self-interaction of the elect ron^&iil 
Correspondingly, the U pcA term in Eq. is known to 
be "made of" two parts: — U po1 = —U p ° l + f/?° , where 
U p ° l represents the magnitude of the potential energy of 
the electron in the polarization field and U po1 the en- 
ergy required to create this polarization ("deformation 
energy"); with the optimal polarization, U%° 1 = U p ° /2. 

Both energy functionals (JIJ and (@J assume that the 
electron and hole energies are measured from the band 
edges. 

It is convenient to factor out dependences on physically 
relevant combinations of parameters by introducing ap- 
propriate units of energy and length. We will choose such 
units based on combinations for the exciton Bohr radius 
and binding energy (effective Rydberg) in 3d: 



ao = 2e?i 2 /me 2 , Ry = e 2 /2eao- 



(5) 



For the exciton problem, Eq. |JT)|. one uses e = too in 
Eq. (J5J and for the polaron problem, Eq. 10} , it is e = e*. 



We use superscript indices "exc" and "pol" to distin- 
guish between the corresponding units JSJ. With all the 
coordinates (x, y for the cylinder) measured in appropri- 
ate ao, one arrives at dimensionless energy functionals: 



cexe 
c 



and 



£ CXC /Ry cxc and £ po1 = £P ol /Ry po1 , where 



£ exc M= dr |V^(r)| 2 -2 / dr 



A)(r) 



(6) 



£ poi m 



dr \Vip(r) 



dri dr 2 



l^(ri)| 2 |V>(r 2 )| 2 

£>o(n - r 2) 



(7) 

The dimensionless -Do( r ) in Eqs. © and Q features the 
same behavior as Eq. J2J but with R replaced by the 
corresponding R$ = R/ag- 



The global minima of Eqs. JJjJ and Q: 



—E po1 , respectively, would determine the dimensionless 
binding energies. As units JSJ already establish the scal- 
ing rules, in what follows we will be comparing _Bg xc and 
E po1 at the same values of Rq. The ratio E po1 /Eq xc would 
have a direct physical meaning of maximum achievable 
when e s 3> Coo and e* ~ eoo in Eq. J2J). As we mentioned 
earlier, this can be a typical situation for many polar 
solvents. 

Before proceeding with the analysis for a cylindrical 
surface, we recall in more detail benchmarks known for 
isotropic d-dimensional systems (Do(r) = |r| = r) briefly 
described in the Introduction. The exact isotropic exci- 
tonic ground state tp(r) corresponding to Eq. |J§} is given 
by the solution of the Schrodinger equation 



-E^ xc and 



TTTCXC 

~ E 



tl> = - 



d 2 j> 

dr 2 



yielding well-known 



ijj(r) oc exp 



2r 
d- 1 



(d — I) difj 2 , 

« i>, 

r or r 



7T1CXC 

b 



(d-iy 



(8) 



With our choice of units, Eq xc = 1 in 3d. 

The polaronic ground-sate wave functions and ener- 
gies corresponding to Eq. J7J) are known from variational 
calculations for d-dimensional systems. It is customary 
in the polaronic literature to express energies in terms 

1/2 

of the coupling constant a c = (me 4 /2e* 2 7i 3 w) and 
phonon frequency u. Note that the combination a 2 hu> 
appearing in results for strong-coupling (adiabatic) po- 
larons corresponds to 2Ry po1 as defined in Eq. (JSJ). Thus 
well-known Pekar's result for the 3d polaron translates 
0.218 and into the same magnitude of the 
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into E ( 

ratio E pol /E$ xc in 3d. 

The 2d polaron has been studied in great detail in 
Ref. E with the best result of E po1 ~ 0.809 for the adia- 
batic case achieved with Pekar-type trial wave functions. 
Since £§ xc = 4 in 2d, the ratio E pol /E$ xc ~ 0.202, only 
slightly smaller than in 3d. The 2d case represents the 
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FIG. 1: The spatial distribution of unnormalized density 
|-!/>(r)| 2 for the exciton (dashed line) and the polaron (solid 
line) in 2d. The polaron wave function is as per Eq. ©. The 
inset shows the corresponding spatial distribution U(r) of the 
potential energy in U = f U(r)dr. 



limit R — > oo for a cylindrical surface and is particu- 
larly important for our analysis. Wc have looked at sim- 
pler one-parametric trial wave functions that would make 
a good representation of the 2d adiabatic polaron and 
found that the wave function 



ip(r) oc 



1 



cosh (ar) 



(9) 



with a ~ 1.674 yields a very good optimization for the 
energy: E^° l ~ 0.804, quite accurate for our purposes. 
Figure n compares the spatial structure of the 2d exciton 
and polaron. 

As both exciton and polaron binding energies diverge 
in pure Id, the following comparison might be of a du- 
bious nature but is still interesting. Specifically, Ref. |6| 
discussed the calculation of the Id polaron in terms of the 
renormalized coupling constant a' c = a c /(d — 1), where 
d — ► 1. The best variational result achieved for the adia- 
batic polaron was E^° l ~ 0.333 {a'^Tiuj). As the diverging 
^-dependence in this expression is the same (d — 1)~ 2 as 
in the exciton case ||SJ|, the ratio of the binding energies 
in Id could then be interpreted as Eq o1 /Eq xc ~ 0.167. 



III. VARIATIONAL ANALYSIS FOR A 
CYLINDER 

When on a cylindrical surface, both exciton and po- 
laron ground states need to be determined numerically. 
Our variational analysis of the energy functionals © and 
has been performed on the following classes of the 
trial wave functions^ For the exciton problem we have 
used three-parametric (a, (3 and 7) wave functions 



ip(x, y) oc exp — (a 2 a; 2 + f3 2 yl + 7 2 ) 



2\V2 



(10) 



The polaron problem is much more computation time 
demanding and we chose two-parametric (a and (3) wave 



ip(x,y) oc 



cosh (cP-x 1 + f3 2 y 2 ) 



2^/2 ' 



(11) 



The functional dependences in Eqs. (|10|l and are 
such that they can recover, in the limit of i?o — * 00, wave 
functions JHJ| and © found for the 2d systems - similarly 
to the earlier exciton calculationsiii*^ 

We explored two choices for the effective coordinate y\ 
in Eqs. ifTU J) and iJTT j) : "arc-based" (as in Ref. EH 

vi = v, (12) 

—ttRq < y < ttRo, and "chord-based" (as in Ref. Hi]) 

2/1 = 2R sin {y/2R Q ) . (13) 

Both choices can be thought of as respectively n>l and 
n = 1 limits of more general 
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ttRq { — sin 

7T 



u 



ttRq 



l/n 



that could be used in future refinements as being more 
flexible in terms of the shape of the wave function pe- 
riodic in the circumferential direction. In this paper we 
resorted to just choosing the best results among obtained 
with Eqs. CE1) and JT3J). 

Wave functions (|10(l and (jl T|) feature two parameters 
a and (3 having the meaning of inverse lengths thereby 
explicitly allowing for anisotropy of the wave function ex- 
tent in the axial and circumferential directions*!* 2 What 
we will later be referring to as quasi-ld results corre- 
sponds to (3 = when the wave functions are uniform 
around the cylinder circumference. 

Main quantitative results of this paper are displayed in 
Figure [3 showing the optimized variational outputs for 
the exciton and polaron binding energies as well as their 



ratio E. 
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Our exciton data is very close to results 
of Refs. [llj and where the reader can find extensive 
discussions. 

Just as in the exciton case, the polaron binding energy 
exhibits very little change from its 2d value (~ 0.8) due 
to the curvature up to Rq ~ 1 (see the inset in Fig. Ufa)) 
where it starts rising, relatively earlier and more rapidly 
than the exciton dependence. This immediately trans- 
lates into a substantial increase of the ratio #p o1 /#o xc 
with decreasing Rq as shown in Fig. H£b). The max- 
imum of the ratio ~ 0.35 is achieved in the region of 
1/Ro ~ 10, after which the ratio starts slowly decreasing 
with 1/Ro- While, of course, larger ratios Eq o1 /Eq* c lead 
to larger relative reductions of the effective activation en- 
ergy Eq xc — 2Eq° 1 , it is quite interesting that the absolute 
value of this activation energy, also shown in Fig. [HTb), 
exhibits a non-monotonic dependence on I/Rq. A region 
around the minimum of this curve indicates specific tube 
sizes where the activation energy would be at its lowest. 

Our variational results have shown that the the depen- 
dence of Eq° 1 on 1/Ro collapses onto the corresponding 
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FIG. 2: (a) The dimensionless exciton and polaron binding 
energies as functions of the dimensionless inverse cylinder ra- 
dius 1/Ro- The inset shows the behavior of the polaron bind- 
ing energy in more detail at small 1/Ro- (b) The ratio of the 
binding energies Eq° 1 /Ef* c of the excitations on a cylinder is 
shown as a thick solid line. A thin solid line displays this ratio 
as it is obtained in the quasi- Id calculation with the effective 
tube potential l|15|l . A thin short-dash line shows the quasi- 
ld result for this ratio if the effective interaction, instead, 
was that of a quantum wire, Eq. 11161 . A thick long-dash line 
displays the dimensionless activation energy _Eg xc — 2Ef° l . 



possible improvements would be inconsequential for our 
main observations of a non-monotonic dependence of the 
ratio E$ ol /E^ c on 1/R and of the magnitude of the 
ratio reaching values well above the 2d value of ~ 0.2. 
We found that a rise of the ratio above 0.3 is obtainable 
even if the polaron wave functions are not specifically 
optimized for a range of given Rq but the 2d optimal 
values of a = [3 ~ 1.67 are used. For very small tubes 
with 1 /itVi > 20, results of exact numerical calculations 
in Ref. 10 complement the picture. 

The transitional region of Rq ~ 1 is likely to be es- 
pecially sensitive to the choice of trial wave functions. 
So Ref. reported a few-percent improvement for the 
exciton binding energy in the region of 1 < l/-Ro ^5 2.5 
with certain trial functions. Similar improvements could 
perhaps be found for the polaron binding energies. Fig- 
ure |3| illustrates the behavior of the polaron functional 
Q in the transitional region, at Rq = 0.65, as a func- 
tion of the variational parameters a and (3. A curious 
feature of the "landscape" in Fig. E^a) is a clear coexis- 
tence of two minima, one corresponding to a polaronic 
state that is uniformly distributed (delocalized) around 
the cylinder circumference, and the other where the cir- 
cumferential distribution is non-uniform. While this ap- 
pears as an interesting possibility, it could as well be an 
artefact of a restricted nature of a specific class of the 
trial wave functions, compare, e.g., visually to the land- 
scape of Fig. |3f b) . A much more careful study would be 
needed to explore the possibility of coexisting polaronic 
states requiring an analysis of actual adiabatic potential 
surfaces. An example of an analogous analysis can be 
found in Ref. l22l where we proved a coexistence of dif- 
ferent polarons in certain quasi-ld systems. 



IV. DISCUSSION 



quasi-ld curve practically right away after the onset of a 
substantial rise in the binding. This is different from the 
exciton case where deviations from the quasi-ld behavior 
persist all the way into the region of 1/Ro > 20. In other 
words, the polaron spreads uniformly around the cylinder 
at much smaller curvatures than the exciton does. We 
note that the quasi- Id variational results for the polaron 
binding using trial wave functions and l|10[l differ 
very little and are in good agreement with our analysis 
in Ref. ^3 where no assumptions have been made about 
the wave function shape and the nonlinear optimizing 
equation has been solved numerically. The thin solid line 
in Fig.^b) shows the ratio Eq° 1 /Eq* c using quasi- Id re- 
sults for the binding energies and its deviation from the 
variational result for a cylinder is entirely due to under- 
estimation of the exciton binding. 

As with all variational calculations, we, of course, can- 
not exclude that some details in the results shown in 
Fig.[21may undergo slight modifications upon further im- 
provements of variational wave functions. Importantly, 



Both in this paper and in Ref. llfj we have shown that 
the strong polaronic effect occurring in sluggish polar en- 
vironments substantially affects the relative energetics of 
closely-bound and well-separated electron-hole pairs on 
a cylindrical surface of nanotubes. The binding energy 
Efr° l of an individual polaron can, in principle, reach as 
much as about 0.35 of the binding energy E^ c of the 
exciton. This would translate into a reduction of the ac- 
tivation energy E^ xc — 2E^° l for exciton dissociation by a 
factor of about 3 from the value E£ KC it would have in a 
non-polar environment with the same value of the high- 
frequency dielectric constant e,^. Note that we have not 
found any additional energy barriers between the exci- 
ton and distant polaron-pair states^ One should expect 
that enhanced separation of charges and a correspond- 
ing luminescence quenching would then result, e.g., in 
experiments probing electric field effects on the lumines- 
cence. Needless to say that additional factors can make 
the reduction magnitude we discussed smaller, see, e.g., 
a comparison in Ref. ^3 of cases with different "electro- 
static conditions". 
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FIG. 3: The behavior of the functional Q at R — 0.65 with 
the trial wave functions il l H as a function of the variational 
parameters a and (3 for two choices of the effective coordinate 
j/i : (a) as in Eq. ifT^jl and (b) as in Eq. 



Polar liquid environments such as many common sol- 
vents may be good candidates to provide conditions nec- 
essary for the strong polaronic effect as they can exhibit 
both high values of the static dielectric constant e s and a 
relatively slow response of the orientational polarization 
(longitudinal relaxation times can be on the order of 1 ps 
and longer) 10,2 «. Quite befittingly, there is an intense on- 
going research effort on various applications of nanotubes 
in contact with such environments and we hope that di- 
rect experimental verifications of our conclusions would 
be possible. We note that a comprehensive mapping 
of luminescence versus absorption spectra of individual 
SWNTs was in fact achieved in aqueous suspensionsi 23 i 24 
Interestingly, numerical estimates in Ref.[l2j indicate that 
the dimensionless radius Rq for a range of SWNTs may 
be close to 0.1, which corresponds to the region of maxi- 
mum E% ol /E^ c ratios in Fig.^b). 

We have demonstrated that the ratio of the binding 
energies E$ ol /E$ xc has a non-monotonic dependence on 
the cylinder curvature 1/R. As argued in the Introduc- 
tion, our particular observation for the cylindrical geom- 
etry can be a manifestation of a more general crossover 
effect that would be common for other structures in 3d 
polar media when the increasing confinement of the elec- 
tron motion causes a "transition" between d-dimensional 
systems, such as between 3d and 2d (quantum wells) or 



between 3d and Id (quantum wires). Basically, the origin 
of this effect can be related to the fact that the spatial 
extent of an unconfmed polaron is larger than the size of 
an unconfmed exciton (see Fig. Q] for the 2d case) thereby 
making the polaron "respond" to the initially increasing 
confinement in a more pronounced way than the exci- 
ton. Only after the exciton experiences a fuller effect of 
the confinement, the ratio of the binding energies starts 
decreasing. Elaborating more on this idea, the inset in 
Fig. n shows the spatial distribution of the potential en- 
ergy terms U: £ — K — U , for unconfmed excitations. 
It is evident that the "longer-range" contributions to the 
polaron potential energy are relatively more important 
than for the exciton. That is why the effect of the in- 
creasing confinement on the polaron is initially relatively 
stronger. In the particular case of the cylindrical ge- 
ometry, the curvature changes remote physical distances 
(chords instead of arcs, Eq. ©) more than it does close 
distances. 

A meaningful parallel can be drawn with the behavior 
of the ratio U/K of the potential and kinetic energy terms 
in confined systems. Virial theorem for the Coulombi- 
cally bound states in unconfmed d-dimensional systems 
states that U/K — 2 independently of d, which, of course, 
also follows directly from scaling of both functionals (JTJ 
and 10} provided that D(r) = |r| = r. Particularly, po- 
larons in such systems are knowni&SSi to satisfy the fol- 
lowing ratios for various energy terms: E po1 : K p ° l : 
U po1 : Ulf : C/p°' = 1:1:2:4:2. Some of these rela- 
tionships are violated in confined systems. As studied in 
Ref.|26j, the virial theorem ratio u cxc /x°^ c for excitons in 
quantum wells and quantum wires is larger than 2 and, 
in fact, a non-monotonic dependence of U e * c /K exc has 
been demonstrated for quantum wells transitioning be- 
tween 3d and 2d limits. We have found a non-monotonic 
behavior of the ratio U/K for both polarons and exci- 
tons on a cylinder as a function of the curvature 1/R. 
In agreement with our qualitative arguments, at smaller 
1/R, this ratio grows much faster for the polaron than 
for the exciton. At large curvatures, however, the trend 
is reversed and the exciton has larger ratios U/K than 
the polaron. 

It is useful to continue a qualitative reasoning by dis- 
cussing the quasi-ld limit of our results, that is, the 
case of stronger but still finite degrees of confinement. 
One can then use the notion of the effective Coulomb 
potentialSfi*2& here as a function of the Id (axial) dis- 
tance x. (The effective Coulomb potentials for 2d can 
be similarly introduced^) Figure 21 compares three po- 
tentials in units such that at very large distances the 
potentials behave as 

V eS (x) = Ro/x. (14) 

In this limit the electron wave function is delocalized 
around the tube circumference and the effective inter- 
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FIG. 4: Distance dependence of the effective Id Coulomb po- 
tentials for charges on a cylinder, Eq. (1151 . solid line, and in 
a quantum wire, Eq. I|16|l . long-dash line, in comparison with 
the original Coulomb interaction 114L short-dash line. 



action becomes that of rings of charge given by 



V eS ( X ) 



7r[(x/i? ) 2 +4] 



where K(m) = J™ (l — to sin 



1/2 



K 



(x/i?o) 2 +4 

-1/2 



d9 is a complete 



elliptic integral of the first kind. If, instead of a tube, we 
dealt with a quantum wire, the electron wave function 
would be delocalized throughout the cross-section of the 
wire with the effective interaction approximated asi 



V eS (x) 



[x/Ro) + 0.3' 



(16) 



Our qualitative arguments, while confirmed by spe- 
cific calculations, do not appear to be restricted to this 
specific situation. We therefore believe that findings of 
larger magnitudes of the ratio E^ ol /E^ c and of a non- 
monotonic dependence of this ratio on the degree of con- 
finement represent a crossover effect that can be com- 
mon to semiconductor nanostructures in 3d polar envi- 
ronments. Further calculations with different structures 
are needed to validate this conjecture and evaluate its 
quantitative aspects, including in non-uniform polar en- 
vironments. 

Another important venue for future research is an as- 
sessment of the activation energy for exciton dissocia- 
tion in confined semiconductors with unequal electron 
and hole masses, m e ^ m^. In this case electron- and 
hole-polarons have different binding energies and, in ad- 
dition, the exciton itself can cause an adiabatic polar- 
ization of the environment A2i The corresponding "pola- 
ronic" corrections to the exciton binding in quantum- well 
wires have, e.g., been studied in Ref. 27. From a general 
standpoint, one would also like to extend the analysis of 
the effects of polar media on the dissociation of confined 
(15) excitons to the intermediate-coupling 2 *^ case. 
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Both effective tube (|15ll and wire (|16fl Id potentials fea- 
ture modifications of the shorter-range interaction from 
the original Coulomb i|F4"|) due to the transverse spread 
of wave functions thereby eliminating a pure Id diver- 
gence of the ground states for excitons and polarons. As 
we discussed above, the role of the longer-distance in- 
teractions is more important for the polaron than for 
the exciton. Since the relative modification of the origi- 
nal Coulomb to the effective potentials is increasing to- 
wards shorter distances, it is then clear that the ratio 
E^ ol /E^ xc in systems with modified interactions should 
be larger than values around 0.2 in systems with the pure 
Coulomb interaction. Moreover, following the same logic, 
one should expect that the larger modification is from the 
original Coulomb distance dependence the larger would 
be Efr° l /E£ xc ratios. Figure^Jshows that deviations from 
the Coulomb dependence for the effective wire potential 
are smaller than for the tube potential (over a relevant 
spatial range). We have performed a quasi-ld variational 
optimization of the polaron and exciton binding ener- 
gies with the tube potential l|15fl . The resulting ratios 
Eq o1 /Eq xc are shown in Fig.^b) with a short-dash thin 
line and indeed smaller than the ratios calculated with 
the tube potential l|15|) . solid thin line in that figure. 
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